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The phase diagram of the first layer of ^He adsorbed on a single graphene sheet has been calculated 
by a series of diffusion Monte Carlo calculations including corrugation effects. As the number of 
C-He interactions is reduced with respect to graphite the binding energy of ''He atoms to graphene is 
approximately 13.4 K/He atom smaller. Our results indicate that the phase diagram is qualitatively 
similar to that of helium on top of graphite. A two-dimensional liquid film on graphene is predicted 
to be metastable with respect to the commensurate solid but the difference in energy between both 
phases is very small, opening the possibility of such a liquid film to be experimentally observed. 



Graphite is a a well known form of carbon, made of two 
dimensional carbon layers glued to each other by interac- 
tions of dispersion type in the z direction, and separated 
by a distance of 3.35 A. Within each of those two dimen- 
sional layers the carbon atoms are located in the nodes of 
a honeycomb lattice, each of them being bound to three 
others by covalent interactions. Even though it is well 
known that to exfoliate graphite is relatively easy, it was 
only recently reported the isolation of a single and stable 
two-dimensional sheet of carbon by mechanical cleavage 
[i,, 2\ ■ This structure is termed graphene and it has been 
predicted to be unstable since the thermal fluctuations 
would make the crystal structure collapse However, 
the experiments show that at least it is kinetically stable 
This novel singular material has already attracted 
the attention of the scientific community, basically for 
its novel electrical properties 4, 5, 6, 7]. 

In this work, we are interested in graphene as a new 
adsorber. Since graphite is a set of graphene layers, one 
would expect a difference in the binding energy of the 
adsorbed species that could lead to a change in their 
phase diagram. We have carried out this analysis for 
^He in the limit of zero temperature, and studying the 
differences between graphene and graphite. To this end, 
we have performed Diffusion Monte Carlo simulations of 
a system of ''He atoms adsorbed on top of a variable 
number of graphene layers, ranging from one to eight. 
This last number was found to be an acceptable model for 
graphite, since all the properties calculated were similar 
for eight and nine layers within the error bars obtained 
from the simulation data. The layers were supposed to be 
parallel to each other, separated by the typical graphite 
distance, and stacked in the A-B-A-B way characteristic 
of this compound. The helium densities were kept within 
the limits of a first layer ( < 0.12 A'^ @). 

We used Diffusion Monte Carlo (DMC) because it is 
able to obtain the true ground state for bosonic systems, 
such as a set of ''He particles Q in the limit of K. 



However, for the technique to work we have to provide 
a reasonable approximation for the ground-state wave 
function, what it is called the trial function, that collects 
all the information known a priori about the system. In 
this work, we used as a first trial wave function 



$(ri,r2, ...,rN) = J|cxp 



i<j 



-'He-He 



(1) 

that depends on the coordinates of the helium atoms 
ri, r2, rjv, and where the first term is the usual Jastrow 
function depending on the helium interatomic distances 
, with 6hc-Ho = 3.07 A^]. For the second part of the 
wave function (\E'(z)), we followed Withlock and collabo- 
rators [l^ , and solved the one-dimensional Schroedinger 
equation describing a single helium atom moving along 
the axis perpendicular to the graphene layers (z) for an 
averaged C-He potential that neglected corrugation. The 
one-body ground state wave function obtained, ^'(z), is 
displayed in Fig. 1. The He-He potential was taken from 
Ref. [llj while the individual C-He interactions were as- 
sumed to be of Lennard- Jones type 1^. This means 



that in our many-body calculations the effects of carbon 
corrugation in the C-He interaction are fully considered. 

Obviously, the former trial function ([T]) is an adequate 
representation for a system with translational invariance, 
i.e., a liquid. To do simulations for a solid, we multiply 
the previous <I'(ri,r2, ...,r]\j) ([1]) by a term of the form 



]^exp-a[(a;i - 



(2) 



where a^sitoi 2/sitc are the coordinates of the crystallograph- 
ical positions around which the '*Hc atoms are localized 
and a is a constant variationally optimized. These po- 
sitions were different for each of the solid phases con- 
sidered: a commensurate -x/S x of surface density 
0.0636 A~^ 141, a commensurate structure (7/16) re- 
ported by Corboz and coworkers [l3| {p = 0.0835 A~^), 
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FIG. 1: Solution for the one-body Schroedinger equation de- 
scribing a single *He atom under an averaged C-He interaction 
in the x,y plane, z indicates distance to such plane, located 
at z = 0. 
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FIG. 2: Energy per atom versus the density for a liquid phase 
of helium atoms on top of one or several graphene sheets. 
Open squares, n — 1; full squares, n — 2; open circles, n — 3; 
full circles, n—4; open triangles, n — S. 



and three commensurate structures suggested to appear 
in the experimental phase diagram of graphite in Ref. 
M (2/5, p = 0.0763 A-2; 3/7, p = 0.0818 A^^ . 3^^/75 
p = 0.0763 A~^) . We also considered several triangular 
incommensurate solids of different densities, obtained by 
varying the helium-helium distance in the x,y plane on 
the (outer) graphene sheet. The optimal values of the 
parameter a in Eq. ^ are a = 0.31 A~^ for all the 
commensurate structures, and ranges from a = 0.15 to 
0.77 A~^ in the case of the incommensurate triangular 
structures. The quality of the different trial functions to 
describe the system can be ascertained by looking at the 
energy variances of every single calculation, in all cases 
of the order of 0.01-0.02 K. 

Fig. 2 displays the energy per ^He atom as a function 
of the helium density for the liquid phase and for the 
number of graphene layers (n) considered. The corre- 
sponding densities were obtained by varying the number 
of helium atoms on top of a fixed simulation cell, a rect- 
angle of 34.43 X 34.08 A^. The use of a simulation cell 
with different size or form did not change the results. We 
can see that there is a very small but significant differ- 
ence between the binding energy values for n = 4 and 
n = 8 (-142.37 ± 0.01 K versus -142.69 ± 0.01 K), for 
the density corresponding to the minimum energy, 0.044 
A~^. This equilibrium density is nearly equal to the one 
for purely two-dimensional liquid '''He (0.043 A~^) [l^. 
The energy differences between structures with the same 
densities for eight and nine carbon sheets are below the 
corresponding error bars and therefore n — 8 can be con- 
sidered the graphite limit. We also found that the curves 
are similar to each other, up to the point to nearly col- 



lapse in one single function when the differences in the 
infinite dilution limit are corrected. The energies of the 
liquid phase at the equilibrium density are shown in Ta- 
ble I. The infinite-dilution energies were obtained from 
fittings to a third-degree polynomials in the density range 
p < 0.02 A-2. The binding energy of a single atom for 
n = 8 is fully compatible with the experimental data 
17 . [l^ . and slightly different from the Path Integral 
Monte Carlo calculations of Ref. 19 (-143.09 ± 0.27 K 
versus -141.64 ± 0.01 K for the present work), probably 
because of the different C-He interactions used. 

In the scenario described above the liquid is not the 
most stable phase at T = K. The energies of the liquid 
phase at equilibrium for different values of n are given in 
Table I, and can be compared with the binding energy for 
the -s/S X commensurate structure. The simulation of 
the X \/3 solid phase has been performed using 120 
atoms in a 44.27 x 42.60 A^ cell. In the Table, we can see 
that in all cases, commensurate solids are more bound, 
but admittedly for very small margins. The ground state 
of ^He on top of any graphene compound is the x VS 
registered phase. This result is stable within plausible 
uncertainties (~ 5%) of the C-He interaction energy scale 
(e): we have verified that the liquid would appear as the 
stable phase by decreasing much more the energy e (20 
%). Our description of the phase diagram agrees with the 
Path Integral Monte Carlo one of Pierce and Manousakis 
[lil [20j for graphite. However, the small difference be- 
tween the energies of the liquid and commensurate phases 
makes also plausible the scenario given by Greywal and 
Busch [3], with a liquid phase of density '-^ 0.04 A~^ as a 
very close metastable state. In the densities between two 
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TABLE I: Energy per atom, in K, for several helium arrangements, n indicates the number of graphene layers considered. See 
further explanation in the text. 



n 


Infinite dilution 


Liquid 




Incommensurate 


1 


-128.26 ± 0.04 


-129.221 ± 0.009 


-129.282 ± 0.007 


-126.6 ± 0.2 


2 


-139.02 ± 0.01 


-139.96 ± 0.01 


-140.067 ± 0.009 


-137.3 ± 0.2 


4 


-141.24 ± 0.09 


-142.37 ± 0.01 


-142.45 ± 0.01 


-139.7 ± 0.2 


8 


-141.64 ± 0.03 


-142.69 ± 0.01 


-142.81 ± 0.01 


-140.0 ± 0.2 
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FIG. 3: Energy versus the inverse of the density for helium 
on top of a single graphene layer. Full circles, liquid phase; 
open circles, incommensurate triangular lattice; open square, 
registered \/3 x \/3 phase. Where not displayed, error bars 
are of the same size or smaller than the symbols. The dashed 
line is the result of a third order polynomial fit to the incom- 
mensurate results. 



stable structures in a phase diagram, the system divides 
itself in patches of the coexisting phases. Those patches 
could be puddles of liquid or clusters of the -\/3 x -s/S 
I9I . 20| surrounded by empty space. 



structure 

The experimental phase diagram of the first layer of 
'^He on top of graphite indicates that at high enough 
density, the stable phase is an incommensurate triangu- 
lar phase In between the y/Sx y/3 registered phase 
and the incommensurate one, different phases have been 
sugg ested 0,[il,[i3, from a domam wall phase (DWP) 
to several commensurate structures of different densities. 
We performed DMC calculations to check the stability of 
these commensurate phases on top of graphene versus a 
triangular arrangement of the same density. The results 
are displayed in Table II. The DWP was not checked be- 
cause in a DMC calculation the atoms in a wall have to be 
located in definite positions. Those positions have to be 
arbitrary (for instance, crossing the simulation cell diag- 
onally or vertically, or in arrays of three atoms instead of 
two) since no experimental information is available. This 
means that the simulation results do not represent the 



DWP phase meaningfully, only giving the corresponding 
energies of a particular atom arrangement including a 
wall. Table II indicates that all the commensurate struc- 
tures have energies per atom larger than the correspond- 
ing triangular phase of the same density. However, and 
except in the case of the 2/5 phase, the differences are 
within two error bars. This means that they could be 
metastable states that could be experimentally observed. 

The phase diagram of ^He on top of a single graphene 
layer can be established with the help of Fig. 3. There, 
we show the energy per "^He atom as a function of the 
inverse of the surface density. The error bars of the in- 
commensurate structure are noticeable larger than in the 
other cases, because they are not simply the statistical er- 
rors of the simulated energies, but the result of averaging 
over different positions of the helium crystallographical 
sites on top of the graphene layer. For each point in the 
figure, we performed four different calculations consider- 
ing similar triangular helium lattices, but displaced a lit- 
tle with respect to each other and averaged the energy re- 
sults. Thus, we take into account the incommensurability 
of solid ■*IIe on top of graphene. The only stable commen- 
surate solid, should be in equilibrium with a triangular 
incommensurate structure. The limits of the coexistence 
phase should have to be determined by a double-tangent 
Maxwell construction. Since the commensurate phase is 
defined by a single density, we can only approximate the 
result, by drawing a line from this point to the one in the 
triangular lattice energy curve that would result in the 
smallest pressure. This imperfect solution gives us that 
the V3 X 73 structure is in equilibrium with a triangu- 
lar phase of ~ 0.08 A~^, in agreement with experimental 
data for graphite Q . In between both equilibrium densi- 
ties, there would be a coexistence zone formed by patches 
of both solids J^], forming a DWP. This entire picture 
for the phase diagram in the n = 1 case is common to 
n — 2,4,8. The energies per helium atom of the Vs x \/3 
phase and the triangular one at 0.08 (equilibrium 
density), are given in Table I (Incommensurate). The 
main difference between graphene and graphite is then 
an offset of ~ 13.4 K in the binding energies of "^He of 
the different compounds. 

All the above results for graphene and graphite were 
calculated without taking into account any three body C- 
He interaction. To do so, one can introduce the so-called 
McLachlan interaction, between a carbon substrate rep- 
resented by a semi-infinite slab and a couple of helium 
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TABLE II: Energy per atom, in K, for several commensurate helium structures. The results labeled "incommensurate" indicate 
the energies for a triangular phase of the same density as the one in the previous line. 



Compound 



Graphene 


-125.81 


± 


0.01 


-126.50 


± 


0.02 


-126.07 


± 


0.01 


-125.89 


± 


0.01 


Graphene (incommensurate) 


-127 


± 0.2 


-126.8 


± 


0.2 


-126.3 


± 


0.2 


-126.0 


± 


0.2 


Graphite 


-139.25 


± 


0.01 


-139.96 


± 


0.01 


-139.54 


± 


0.01 


-139.33 


± 


0.01 


Graphite (incommensurate) 


-140.5 


± 


0.2 


-140.2 


± 


0.2 


-139.7 


± 


0.2 


-140.0 


± 


0.2 


Graphite (+ McLachlan) 


-138.75 


± 


0.03 


-139.50 


± 


0.02 


-139.02 


± 


0.01 


-138.81 


± 


0.01 


Graphite (+ McLachlan incommensurate) 


-140.1 


± 


0.2 


-139.6 


± 


0.2 


-139.3 


± 


0.2 


-138.9 


± 


0.2 



atoms on top of it 2l|. This means that the McLachlan 
term only can be applied meaningfully to graphite. The 
model of graphite used in our calculations is not a slab 
but a stack of graphene sheets, implying that the use of 
this term is only an approximation. In previous calcula- 
tions 1^ lOi 12| J this term was found to favor the liquid 
versus the commensurate -s/S x VS structure. To check 
its influence on the energy of the system, we performed 
DMC calculations on the ^/3 x ^/3 and liquid phases on 
top of graphite in the very same conditions given above. 
The binding energy for the commensurate solid was - 
142.44 ± 0.01 K versus -142.81 ± 0.01 K of Table I, i.e. 
a difference of 0.37 K. The density minimum for the liq- 
uid phase changed from 0.044 to 0.041 A"^, with 
a binding energy of -142.49 ± 0.01 K versus the -142.69 
± 0.01 K of Table L This would mean that the liquid 
is the stable phase for '*He on top of graphite, but with 
a difference even smaller than in the above calculations. 
We also repeated the calculations including this term for 
all the registered phases suggested above. The results 
are given in Table II. We found that the McLachlan in- 
teraction change the binding energies but not the fact 
that they are metastable with respect to the incommen- 
surate triangular phase. The results are not applicable 
to graphene, since the model of a thick carbon slab does 
not apply to that substrate. 

Summarizing, we have performed diffusion Monte 
Carlo calculations of ^He adsorbed on graphene for the 
first time using an accurate He-He interatomic potential. 
Our results show that the ground state corresponds to a 
•\/3 X VS commensurate solid. However, the difference in 
energy between a metastable liquid film and the solid is 
very tiny opening the possibility of observing experimen- 
tally a two-dimensional superfluid liquid phase. Other 
commensurate phases have also been studied and found 
to be metastable, both in graphene and graphite. The 
phase diagram of ^He on graphene is qualitatively equal 
to the one of ^He on graphite that we have obtained by 
adding graphene sheets up to ti = 8. 
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